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Received 

Abstract. We present optical and near-infrared (NIR) photometry of the bright afterglow of GRB 020813. Our data span from 
3 hours to 4 days after the GRB event. A rather sharp achromatic break is present in the light curve, 14 hours after the trigger. 
In the framework of jetted fireballs, this break corresponds to a jet half-opening angle of 1.9° ± 0.2°, the smallest value ever 
inferred for a GRB. We discuss our results in the framework of currently available models, and find that they have problems in 
explaining the joint temporal and spectral properties, and in particular the slow decay before the break. 
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1. Introduction 

Since the discovery of the first afterglow of a gamma-ray burst 
(GRB; Costa et al. nWTI van Paradjis et al. lTWTt . our knowl- 
edge of these mysterious explosions has rapidly increased. The 
huge energetics implied by their cosmological distance (e.g. 
Metzger et al. ll997> has severely constrained existing models, 
leading to the standard internal/external shock scenario (e.g. 
Piran J^99, Meszaros 2002 and references therein). According 
to this model, the prompt gamma-ray emission is produced by 
internal collisions in a relativistic blastwave, while the after- 

Send offprint requests to: S. Covino e-mail: covinoOmi . astro . it 

* Based on observations partly made with ESO telescopes at the 
Paranal Observatories under programme Id 69.D-0461 and with the 
Italian TNG telescope under programme TAC 8_01(47). 



glow originates by the slowing down of the same fireball in the 
surrounding medium. 

An important question still remains unsettled: are the fire- 
balls spheres or jets? In the latter case, a steepening in the 
light curve is expected when the bulk Lorentz factor of the fire- 
ball equals the inverse of the opening angle of the jet (Rhoads 
1999 1. This steepening must occur at all frequencies and at 
the same time ('achromatic' break). To date, there are several 
possible examples of such behaviour in the optical band. The 
most convincing cases are perhaps GRB 990510 (e.g. Israel 
et al. Harrison et al. 1999 Stanek et al. [D^^ and 

GRB 010222 (e.g. Masetti et al. 2001 , Sagar et al. l200Tl Stanek 
et al. 2001 1, for which the break times were ~ 1.5 days and 
~ 0.7 days respectively. 
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Fig. 1. The field of GRB 020813 (~ 1.5' X 1.5'), as imaged by 
TNG-DOLORES on Aug 15 UT. The afterglow is marked by 
an arrow. This is a composite image from B, V, and R filters. 
North is up. East is left. 



GRB 0208 13 was detected on 2002 August 13 by the 
HETE-2 spacecraft at 2:44:19 UT (Villasenor et al. l20^ . It 
was a bright event, lasting more than 125 s, with a fluence 
of ~ 3.8 X 10"^ erg cm"^ (25 H- 100 keV) as reported by 
Ulysses (Hurley et al. 2002). The rapid coordinate dissemi- 
nation and the small HETE-2 error box (4' radius) allowed 
Fox, Blake & Price (2002) and Gladders & Hall (.2002a) to 
independently identify a bright optical counterpart with posi- 
tion 0-2000 = 19H6'"4P.9, ^2000 = -19°36'04'.'8, just 1.9 hours 
after the trigger The spectroscopic redshift, determined with 
KECK 1-LRIS, was z = 1.255 (Barth et al HUU^ . Subsequent 
observations by various groups allowed a preliminary sampling 
of the light curve, which exhibited a steepening in the optical 
bands beginning a few hours after the GRB trigger (e.g. Bloom, 
Fox & Hunt 2002 1. Polarimetric observations were carried out 
by Barth et al. (2003 1 and Covino et al. (2002 1, yielding pos- 
itive polarization signals at the level of P - (2.9 + 0.1)% and 
P = (1.17 + 0.16)%, 4.7 hours and 23 hours after the burst, re- 
spectively (both measurements are uncorrected for the small 
Galactic-induced polarization). The afterglow was also ob- 
served in the X-ray band (0.6 -H 6 ke V) by Chandra, as a bright 
fading source with an average flux F ~ 1 .9 x 10"'^ erg cm"^ s"' 
(Butler et al. l2003> . A positive radio detection was obtained by 
Frail & Berger ("2002^, with Fy = 300 /Jy at 8.46 GHz. Only 
upper limits were reported in the millimeter region (Bremer & 
Casti-o-Tirado l2002l Bertoldi et al. l2002t . 

2. Data collection, reduction, and analysis 

We began observing the optical afterglow of GRB 0208 13 
(Fig.0 on Aug 13.99 UT, with the ESQ VLT-UT3 (Melipal), 
equipped with FORS 1 and a V filter in imaging polarimet- 
ric mode; during the subsequent nights (Aug 15, 16, and 17), 



we monitored the optical/NIR light curve using the 3.6 m 
Telescopio Nazionale Galileo (TNG) in the Canary Islands. 
UBVRI and JHK data were acquired with the DOLORES 
and NICS instruments respectively. Additional NIR images 
(J and K) were taken with the AZT24 1.1 m telescope in 
Campo Imperatore (L'Aquila, Italy) on Aug 13. We then added 
to our photometric set the publicly available data' of Mike 
Gladders and Pat Hall (.2002b i. acquired with the Baade 6.5 m 
(Magellan 1) telescope. Standard stars were imaged during the 
night of Aug 16 with the TNG. Data reduction and calibra- 
tion were carried out following standard procedures, as imple- 
mented in the Eclipse package (version 4.2.1; Devillard ll997> . 
while photometry was performed by means of GAIA^. All pho- 
tometric data, shown in Fig.|2l are posted online''. 

We modelled the afterglow light curve with the functional 
form first proposed by Beuermann et al. ( !1999> : 



[(f/fb)**. +(f/fb)'^^2]i/^' 

here Fy is the (monochromatic) flux at frequency v, and t is the 
time elapsed from the GRB; this expression reduces to Fy oc 
r^' for f « fb and Fy oc r''' for f » fb; fb is the break time 
at which the transition between the two regimes occurs, and 
k > describes how fast this transition takes place (the larger 
k the sharper the transition); is a normalization constant. 
We started by fitting the V, R and / bands, which are sampled 
closely enough to constrain fb and 6i. To within the errors, the 
inferred parameters, summarized in Table [0 were determined 
to be the same between the dififerent bands, even if there is 
a marginal hint (not statistically significant) for a correlation 
between fb and the frequency. We caution that 62 could also 
be larger, since our sampling does not extend very much after 
fb. By fitting simultaneously the whole dataset, including U, B 
and NIR data, we got our best values. As shown in Fig. |2] at 
early times the decay was fairly slow, similar to that seen in 
other bursts with an early break (e.g. GRB 010222: Masetti et 
al. l2001> . After fb, the slope increased by 0.65+0.05. In all fits, 
the parameter k, though poorly constrained, was always > 5, 
indicating a very sharp transition. This is confirmed by the fact 
that early data (f < 5 h) are well fitted by a single powerlaw, 
showing that the steepening had not yet started at this time. 
This is in contrast with the finding of Gladders & Hall ( 2002ct . 
who claim a steepening in the /-band between 3 and 5 hours 
after the GRB. 

To further check if the break is really achromatic, we stud- 
ied the spectral properties of the afterglow, using our photo- 
metric multiband data points, both before and after the break 
(Fig- Ell. We first corrected for Galactic reddening, assuming a 
neutral hydrogen column density'* of A^h = 7.52 x 10^" cm"^, 
corresponding to a V-band absorption Ay - 0.42 (Predehl & 
Schmitt 1995^). Individual powerlaw fits to each dataset yield 
the spectral indices reported in Table |2 (Fy oc v "), consis- 
tent with a remaining constant during all the observations. This 

' ftp://ftp.ociw.edu/pub/gladders/GRB/GRB(S)20813 . 
^ http://star-www.dur.ac.uk/~pdraper/gaia/gaia.html . 
' http : //www .merate .mi . astro . it/~malesani/GRB/828813 
http://heasarc.gsfc.nasa.gov/docs/corp/tools.html . 
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Fig. 2. Optical and NIR light curves of GRB 0208 1 3 afterglow, 
with our best fit according to Eq. [2 (dashed lines); refer to 
Table[l]for fit parameters, and to the web page"* for the numer- 
ical values. For the first night (f < 5 h), error bars are smaller 
than symbols and are not plotted. The thick solid line shows the 
(arbitrarily normalized) X-ray decay. 

Table 1. Fitted parameters describing the temporal evolution 
of the afterglow. 6i and 62 are the early- and late-time temporal 
slopes, and fb is the break time (measured in hours since Aug 
13.1141 UT); see Eq.[T] Errors are 1-cr. 



Band 


Si 




62 




tb(h) 




V 


0.84 


±0.11 


1.34 


±0.06 


11.9 ± 


6.4 


R 


0.76 


±0.05 


1.46 


±0.04 


13.7 ± 


1.2 


1 


0.80 


±0.03 


1.43 


±0.12 


16.7 ± 


1.7 


All 


0.78 


±0.02 


1.44 


±0.04 


14.2 ± 


0.8 



Table 2. Time evolution of the spectral index (Fy v""). 
Galactic (but not intrinsic) extinction has been accounted for 
by assuming Ay = 0.42. 



UT 


# 


Telescope 


Bands 


a 


13.28 


A 


Baade 6.5m 


VRl 


1.06 ±0.06 


13.31 


B 


Baade 6.5m 


VRI 


1.00 ±0.05 


13.88 


C 


AZT 1.1m 


RJK 


0.85 ±0.14 


14.04 


D 


Baade 6.5m 


BVRI 


1.15 ±0.05 


15.04 


E 


TNG 3.5m 


UBVRIJHK 


1.01 ±0.04 


16.01 


F 


TNG 3.5m 


BVR 


0.72 ± 0.22 


All 








1.04 ±0.03 



is a confirmation that the break is indeed achromatic. A si- 
multaneous fit to the full dataset (six spectra) gives the best 
value a = 1.04 + 0.03. The colors are B - V = 0.32 + 0.05, 
V-R^ 0.39 + 0.06, R-I ^ 0.45 + 0.10, within the range of 
the sample by Simon et al. (2001 1. The chi-square of the fit is 
= 82 for 17 degrees of freedom. Such large value is due to 
the use of data coming from different telescopes. 
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Fig. 3. a) Photometric spectra of GRB 028013 afterglow at var- 
ious times (see Table|2}, together with the best fit (dotted lines). 
Data are corrected for reddening in the Galaxy but not in the 
host. Error bars are smaller than the symbols and are not plot- 
ted, b) TNG spectrum on Aug 15.04 UT; open points and dot- 
ted line are corrected only for Galactic reddening, while filled 
points and dot-dashed line are corrected also for reddening in 
the host (see text). 

3. Discussion 

The afterglow of GRB 0208 13 represents a new case in 
which the optical light curve shows an achromatic break. The 
break was quite sharp and occurred early, similar to that of 
GRB 010222. The slope after the break (~ 1.4) was typical of 
afterglows at times between 1 and 10 days after the GRB. Since 
observations are often performed during this range of days (and 
only rarely earlier), it is possible that many such breaks have 
been missed in the past. An intriguing example of this situa- 
tion is GRB 010921, whose light curve might have displayed 
two breaks, the first only constrained by an early LOTIS upper 
limit to be at less than ~ 1 day (Park et al. 2002 1, and the sec- 
ond determined by means of HST observation after ~ 35 days 
(Price et al. HUU^ . 

There are various interpretations to explain the presence of 
an achromatic break. If the blastwave producing the GRB is 
collimated, the break is expected to occur when the inverse of 
the bulk Lorentz factor of the ejecta F equals the half-opening 
angle of the jet )?jet (Rhoads 1999i. In this context. Frail et al. 
(1200 li found that the net energy released by GRBs, after cor- 
recting for beaming, is approximately constant for a number of 
events, to within a factor of ~ 2. This result is based on the ob- 
served correlation between break times and measured isotropic 
energies. Using the fluence 5 = 3.8 x 10"^ erg cm"^ reported by 
Hurley et al. (2002 1 in the (25 H- 100) keV band, and applying a 
bolometric correction of 5.5 + 3.5 (calculated in the same way 
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Fig, 4, Radio to X-ray spectral energy distribution of the af- 
terglow of GRB 020813, 30 h after the trigger. Optical (BVRI 
bands) and X-ray data have been computed at the epoch of the 
radio detection by interpolating between different spectra and 
using the appropriate decay law. Open symbols are corrected 
only for Galactic extinction, while the additional host com- 
ponent has been taken into account for filled symbols. Radio 
detection together with the millimeter upper bounds limit the 
spectral slope {Fy oc v"") to be a > -1/3 at GHz frequencies 
(in particular a self-absorbed spectrum is ruled out). In the re- 
gion between the milUmeter and optical bands, the drawing is 
only indicative. 



as in Bloom, Frail & Sari l200lT l. we get an isotropic energy^ of 
(9.8 + 6.2) X 10^^ erg, ranking third in the sample of bursts with 
known redshift after GRB 990123 and GRB 000131. Using the 
formalism described by Frail et al. (2001), the opening angle 
of the jet is )9jet = (1.7° ± 0.2°)nJ^^ri2^ , where n_i is the exter- 
nal number density in units of 10"' cm"-' and tj - 20% 7720 
is the gamma-ray production efficiency. This is the smallest 
value ever reported for a GRB. The beaming-corrected energy 
is hence (5.4 + 2.5) x IQ^^n^J^ii^i^ erg, close to the value found 
by Frail et al. (1200 1> in their sample. 

In a different interpretation, an early break can be caused 
by the onset of the nonrelativistic phase, when F drops to unity 
(e.g. Dai & Lu ll999> . This interpretation has been proposed 
for GRB 010222, which shows spectral and temporal proper- 
ties similar to the case presented here (Masetti et al. 2001). 
In order to stop the fireball so quickly, however, an extremely 
high external density n ~ 10** cm"^ is required (Panaitescu 
& Kumar 2000). Since the radius of the fireball at f = fb is 
R ~ 10'^ cm, the surrounding material must in this case have a 
hydrogen column density of A^h ~ nR ~ 10^"^ cm"^, and hence 
it must be completely ionized in order not to conflict with the 
X-ray observation, which does not show a high value of A^h at 
the redshift of the host (Butler et al. 2003). These parameters 
are typical of a supernova remnant (e.g. Vietri & Stella 1998 ). 
If this interpretation is correct, however, the radio spectrum is 
expected to be heavily self-absorbed; this is not the case, as can 
be seen in Fig.|3 



^ We adopt a standard cosmology: ho = 0.65, fl„, = 0.3, Qa = 0.7. 



Table 3. Predicted slopes and comparison with our observed 
values 5i = 0.78+0.02 and ^2 = 1.44+0.04 for different models 
(all in the framework of a spherical relativistic fireball). Vi and 
Vc are the standard injection and cooling frequencies (defined 
e.g. in Sari et al. 1998J . Vo spans the whole observed band, from 
optical (~ lO'"* Hz) to X-ray (~ 10'** Hz). 



Frequencies 


Ambient 


Electron 


Predicted 


Discrepance 


order 


density 


index p 


slope S 


<5i 


S2 


Vc < Vi < Vo 


uniform 


1.6±0.1 


0.92+0.02 


7.5o- 


13o- 




wind 


1.6+0.1 


0.95+0.01 


9.0cr 


120- 


Vi < Vo < Vc 


uniform 


2.6+0.1 


1.2+0.1 


4.3cr 


2.4a- 




wind 


2.6+0.1 


1.7+0.1 


9.3o- 


2.60- 


Vi < Vc < Vo 


uniform 


1.6+0.1 


0.92+0.02 


7.5o- 


13o- 




wind 


1.6+0.1 


0.95+0.01 


9.O0- 


120- 



To extract further information on the physical parameters 
of the explosion, we attempted to model the broad-band spec- 
tral energy distribution (SED) of the afterglow, including the 
radio detection (Frail & Berger 2002 ), millimeter upper lim- 
its (Bremer & Castro-Tirado l2002l Bertoldi et al. 12003, and 
X-ray data (Butler et al. '2003V The extrapolation of the opti- 
cal spectrum lies well below the value observed in the X-ray 
band, which has moreover a harder spectrum (Fig. ^ dashed 
line). This could indicate that such emission constitutes a dif- 
ferent component (due e.g. to the Compton process, such as in 
GRB 000926; Harrison et al. (20011 ). Alternatively, additional 
extinction could be responsible for the optical/X-ray mis- 
match; we therefore fitted again all optical/NIR datasets with 
a powerlaw allowing for additional absorption at z = 1 .254, 
adopting the extinction law of Cardelli, Clayton & Mathis 
( 1989). The best fit yields a = 0.85 + 0.07 with Av(host) = 
0.12 + 0.04, a small but non-negligible amount. We gstx^ -51 
for 14 degrees of freedom. The statistical improvement (with 
the respect to the case of no additional extinction) is admittedly 
low; however, in this case the optical spectrum becomes harder, 
and a single component can account for both the NIR/optical 
and X-ray emissions (see Fig. |4] shaded region). Moreover, 
the spectral indices in the two bands are strikingly similar: 
ax = 0.85 ± 0.04, Copt = 0.85 ± 0.07. Last, we note that 
(independently on any extinction) the decay slopes are very 
well constrained, and yet in full agreement: Sx = 1.38 + 0.06, 
^opt = 1 .44+0.04. This seems to be a fine tuning if the two com- 
ponents have a different origin. Using Chandra data, Butler et 
al. (l2003t have reported no excess column density with respect 
to the Galactic value. However, our Ay corresponds to a mod- 
est A^H = 1.95 X 10^" cm"^ (assuming a Galactic gas to dust 
ratio), a very small column to detect at a redshift of z = 1 .254 
(the observed value is reduced by a factor ~ (1 + z)^ ~ 10). In 
the following, we assume that indeed NIR/optical and X-ray 
emissions constitute a single component. 

Since ax < 1, the spectrum is hard, with its peak fre- 
quency (in vFy) lying above the X-ray band. If the emission 
we see is synchrotron by a powerlaw distribution of electrons 
(N(y) oc j^P), there are two possibilities, depending on whether 
the emitting particles are cooling rapidly or not. In the first case, 
p - 2a - 1.6 + 0.1 ('flat' distribution), while in the second 
case /:> = 2a; + 1 - 2.6 ± 0.1 (a more conventional value). 
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Theoretical models predict several relations between the tem- 
poral and spectral slopes, depending on a number of assump- 
tions such as the external density profile, cooling regime, and 
dynamical conditions. Before the achromatic break, however, 
the fireball should always follow a relativistic spherically sym- 
metric evolution, and the predictions are robust and easy to 
check (e.g. Panaitescu & Kumar l2000l Dai & Cheng '^OOTi. 
The only delicate issue regards the flat distribution case (when 
p < 2), where the evolution of the spectrum is sensitive to the 
high-energy cutoff of the electron distribution, the behaviour of 
which is not well understood; we adopt the prescription given 
e.g. by Moderski et al. (110001 and Dai & Cheng ( EUOn . Table|3 
shows all possible cases, presenting different ambient media 
and spectral shapes. None of them can account for the low 
value 6] = 0.78 + 0.02. To eliminate such mismatch, a pos- 
sible explanation is to invoke a strong Compton cooling (ei- 
ther on self-synchrotron or external photons), whose emission 
is confined at high energies (> 10 keV). In this case, the cool- 
ing frequency should decrease slower with time, or even in- 
crease, and so the decay would be flatter above this frequency. 
Alternatively, a refreshed shock model (Sari & Meszaros l2000t 
could also alleviate the discrepance, as recently proposed by 
Bjornsson et al. (2002) for GRB 010222. 
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